AD-A035  050 


UNCLASSIFIED 


AO 

A035050 


ARMY  MISSILE  RESEARCH  DEVELOPMENT  AND  ENGINEERING  LAB--ETC  F/G  17/9 
QUASI -COHERENT  DETECTOR  PERFORMANCE  WITH  APPLICATION  TO  CW  RADA— ETC (U) 
OCT  76  T H GEE*  R R BOOTHt 

RE-77-1  NL 


TECHNICAL  REPORT  RE-77-1 


QUASI -COHERENT  DETECTOR  PERFORMANCE  WITH 
APPLICATION  TO  CW  RADAR 


Thomas  H.  Gee 
Virginia  Military  Institute 
Lexington,  Virginia 


Robert  R.  Boothe 
Advanced  Sensors  Directorate 

US  Army  Missile  Reaaarch,  Development  and  Engineering  Laboratory 
US  Army  Missile  Command 
Redstone  Arsenal,  Alabama  38808 


t - 

♦ 


[f . 

I i1 

f <•  t 

r 

K 

■ 


L 


UNCIASSTFTKn 


SECURITY  CLASSIFICATION  OF  THIS  FACE  (Whtn  Dote  Entered) 


REPORT  DOCUMENTATION  PAGE 


|l.  REPORT  HUMBER 

2.  GOVT  ACCESSION  NO. 

in 

ft} 

4.  TITLE  C«nd  Submit) 


'J 


'Quasi -Coherent  Detector  Performance  with  Application' 
to  CW  Radar,  


7.  AUTHOR!*) 


Thomas  H . /^eefvirginia  Military  Institute 
Robert  R. /Boothe)  Adv.  Sensors  Dir. , MICOM 

•TTERFOmSHf^ROANizATION  NAME  AND  AOORESS 

Commander 

US  Army  Missile  Command 
Attn:  DRSMI-RE 

Redstone  Arsenal,  Alabama  35809 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 
1.  RECIPIENT’S  CATALOO  NUMBER 


s./type  of  report  a period  covereo 


rechnlcal 


PERFOI 

RE-77-1 


mjkpgt  ttJ 


NUMBER 


4.  CONTRACT  OR  GRANT  NUMBER!*) 

632303J 11. 21401 


AMCMS  Code  6323 

6^1W36jyS2fT 


Element  ^ede  m a 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  • WORK  UNIT  NUMBERS 


tl.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

Commander 

US  Army  Missile  Command 

Attn:  DRSMl-RPR 

Redstone  Arsenal,  Alabama  35809 


X — •J-U*-StansT  niT»Hr 

( yyji  dctmmmiG  / 
\tlr\  ii  iiumwiap>>AGEr 

50 


14.  MONITORING  AOENCY  NAME  4 AOORESS!!!  rffffarant  from  Controlling  Office.) 


IS.  SECURITY  CLASS,  !«f  0,1  e raport) 

Unclassified 


IS*.  OECLASSI FIC  ATI  ON/ DOWN  GRADING 
SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  (ml  0,1m  Raport) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ml  Oim  mbmtrmet  entered  In  Block  30,  If  dlHmrmnl  from  Raport) 


IS.  SUPPLEMENTARY  NOTES 


IS.  KEY  WOROS  (Continue  on  rereteo  elde  II  nmemmmmtr  mod  Identify  h y Mock  number) 

Coherent  processor 
CW  radar 

Doppler  — phase  processing 


30.  ABSTRACT 


I Identify  hr  block  , 


JThe  evolution  of  a qua  si -coherent  detector  having  application  to  CW  radar  is  pre- 
sented. The  detector  Is  configured  as  parallel  correlator  channels  which  use  digital 
techniques  to  process  video  signals.  Doppler  frequency  channels  and  channels  repre- 
senting the  Initial  phase  of  the  received  signal,  or  return,  are  employed  to  approximate 
a coherent  detector.  A distinguishing  feature  of  the  qua  si -coherent  detector  lies  in  the 


DO  I'jJTn  U73  EDITION  OF  I MOV  ••  IS  OBSOLETE 


UNCLASSIFIED -Jc 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Rm  Dmtm  Unteted'cyf1 

Ifo  3 0 ? ' 


■ ... . kJessf 


> t- 


UNCLASSIFIED 


mas 1 caa  an 


ABSTRACT  (Concluded) 

use  of  a rectangular  pulse  train  as  the  local  signal  which  Is  correlated  with  the  return. 
The  rectangular  pulse  structure  provides  computational  efficiency,  while  being  opti- 
mized with  respect  to  signal -to-nolse  ratio. 

Analytical  results  of  particular  Interest  pertain  to  the  comparative  performance 
of  the  qua  si -coherent  detector  and  the  optimum  coherent  detector,  the  latter  providing 
a detection  performance  benchmark,  hi  the  absence  of  errors  in  estimation  of  doppler 
frequency  and  Initial  phase,  it  is  shown  that  the  detector  performs  about  0. 35  dB 
Inferior  to  optimum.  With  limitations  of  a finite  number  of  processor  channels,  trade- 
offs between  implementation  complexity  and  performance  in  the  presence  of  estimation 
errors  are  discussed.  Criteria  for  design  of  the  qua  si -coherent  detector  are  presented. 


UNCLASSIFIED 


security  cl ASSir«c**» 


Dmtm 


■ 


CONTENTS 


Page 


I.  INTRODUCTION 3 

H.  PROCESSOR  CONFIGURATION 4 

HI.  PROCESSOR  ANALYSIS 6 

IV.  DESIGN  CONSIDERATIONS 22 

V.  SUMMARY  AND  RECOMMENDATIONS 24 

Appendix.  DETECTOR  PERFORMANCE  FOR  VARIOUS  SIGNAL 
PROCESSOR  CONFIGURATIONS  AND  THRESHOLD 
COMPARISON  SCHEMES 27 


ACCESSION  W 

tms 

DOC 

OHMKOUHCEft 

JUSIIHCMIW- 


IT 


Whitt  Stdto 
lull  acuta 


IT 


OmilCUTMU/MMlttlUTT  CHI 

uM-i.-aT-tflHttaLI 


ant. 


□ □ 


I. 


INTRODUCTION 


Radar  returns  are  generally  characterized  as  signals  which  possess 
doppler  dependent  frequency  and  unknown  phase.  Such  signals  can  be  non- 
coherently  processed  to  yield  optimum  detection  performance  within  the  con- 
straints of  unknown  phase.  However,  If  these  constraints  are  removed,  the 
performance  of  such  a system  is  Inferior  to  that  of  a coherent  processor  wherein 
the  phase  Is  exactly  known.  A signal  processor  configuration  is  considered  in 
this  report  which  is  an  approximation  to  the  optimum  coherent  processor  in  that 
it  employs  a number  of  phase  "bins,"  one  of  which  approximates  the  phase  of 
the  return  signal.  In  this  context,  this  signal  processor  might  be  termed 
"qua  si -coherent. " 

The  work  reported  here  was  undertaken  in  support  of  an  on-going  Army 

program  for  development  of  a "quiet"  radar  system.  * The  signal  format  of 
interest  is  a pseudo-noise  (PN)  coded  CW  signal,  hi  particular,  the  analyses 
and  results  which  are  presented  pertain  to  biphase,  PN  coded  signals.  Maxi- 
mal length  PN  code  sequences  are  assumed.  Performance  of  a signal  processor- 
detector  configuration  as  referred  to  in  this  report  is  defined  as  a measure  of 
detection  probability  for  a given  probability  of  false  alarm,  with  signal -to- 
noise  (S/N)  ratio  acting  as  a variable  parameter. 

Accounting  for  range,  doppler  and  phase,  a potentially  large  number  of 
processor  channels  might  be  required  to  achieve  near  optimum  performance. 
Objectives  of  this  reported  work  include  determination  of  the  tradeoffs  which 
are  available  between  implementation  complexity  and  performance.  It  is 
envisioned  that  many  of  the  processor  operations  can  be  best  accomplished  by 
digital  means.  For  this  reason,  attention  is  given  to  configuring  a processor 
which  is  conducive  to  the  use  of  digital  techniques. 

Evolvement  of  the  quasi-coherent  detector  is  outlined  in  the  discussion 
which  follows.  A series  of  modifications  to  the  basic  correlator  configuration 
as  employed  to  process  video  signals  are  described  and  analyzed.  It  is  shown 
that  with  some  restrictions/ modifications  of  the  conventional  correlator  con- 
figuration can  be  effected  to  yield  near  optimum  detector  performance.  Still, 
the  ability  is  preserved  to  employ  digital  processing  techniques  which  are 
simple,  fast,  and  within  the  state-of-the-art. 


Boothe,  R.  et  al. , Quiet  Radar/  Missile  Guidance  Technical  Assessment. 
Final  Report,  US  Army  Missile  Command,  Redstone  Arsenal,  Alabama, 
TR-RE-76-12,  September  1975. 
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II.  PROCESSOR  CONFIGURATION 


The  signal  processor  configuration  of  interest  here  is  shown  in 
Figure  1.  This  figure  illustrates  LxMxN  processor  channels  which  perform 
correlation  operations  between  the  return  signal  and  locally  generated  wave- 
forms. There  are  L range  channels,  where  L is  the  length  of  the  PN  sequence 
employed  to  code  the  transmitted  signal.  This  number  is  determined  in  part  by 
operational  requirements  of  the  radar  system,  hi  addition,  there  exist  M 
doppler  channels  per  range  channel.  Associated  with  each  doppler  channel, 
there  are  N phase  channels.  Given  the  number  L,  it  is  desired  to  determine 
the  number  Mxn  channels  required  to  achieve  a specified  level  of  perform- 
ance for  the  detector.  The  primary  objectives  of  the  analyses  which  follow 
are  to  assess  detector  performance  as  a function  of  the  number,  M x N,  of 
processor  channels  and  S/  N ratio,  and  to  determine  the  interdependence  of 
values  of  M and  N.  It  is  desirable  to  minimize  M x N while  maintaining  an 
acceptable  level  of  performance. 


% 


RANGE 


Figure  1.  Processor  channel  structure. 

Restricting  attention  to  one  particular  doppler-phase  channel,  the  video 
processor  configuration  shown  in  Figure  2 forms  the  basis  for  further  con- 
sideration. With  reference  to  this  figure,  a PN  coded  biphase  signal  with 
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Figure  2.  Processor  block  diagram 


additive  white  gaussian  noise  is  input  to  the  mixers  shown.  The  received 
signal  is 


where  the  noise  i3 


and  c(t)  represents  a maximal  length  PN  code  sequence  with  values  ±1.  Doppler 

frequency  shift  is  denoted  by  o> The  mixer  outputs  are  filtered  to  produce 

d 

in-phase  and  quadrature  (I  and  Q)  channel  video  signals.  Figure  2 indicates  that 
the  PN  code  is  removed  at  the  mixers  by  multiplying  by  a shifted  local  replica 
of  the  code  sequence,  c(t  + r),  where  r takes  on  discrete  values  corresponding 
to  each  range  bin  or  channel,  tt  is  not  necessary  to  decode  at  this  point;  but 
rather,  decoding  could  be  accomplished  at  video,  using  digital  techniques. 


The  video  signals  In  the  I and  Q channels  of  the  processor  now  have  the 


! 

\ 


form 


x(t)  = Ac(t)c(t+  t)  cos  + c(t+  r)n^(t) 


(3) 


and 

y(t)  = Ac(t)c(t+ r)  sin  (w^t+ e)  + c(t+ T)n^(t)  . (4) 

hi  the  in-range  channels,  r equals  zero,  and  ideally 

c(t)c(t+  t)  ■ 1 (5) 

These  signals  are  next  multiplied  by  locally  generated  waveforms,  denoted  by 

g (t)  and  g (t).  The  locally  generated  waveforms  for  each  doppler-phase 
x y 

channel  represent  doppler  frequency  and  phase  estimates  of  the  received  signal, 
hi  this  report,  the  three  multiplier  waveforms  shown  in  Figure  3 are  considered. 
The  sinusoidal  waveform  is  considered  because  it  provides  optimum  perform- 
ance if  the  frequency  and  phase  estimates  are  exact,  while  the  square  wave  and 
rectangular  pulse  waveform  are  conducive  to  the  use  of  digital  signal  process- 
ing techniques. 

The  I and  Q channel  signals  are  summed  as  indicated  by  Figure  2,  and 
the"  input  to  an  integrator  to  complete  the  correlation  process.  Finally,  at 
time  T,  the  integrator  is  sampled  and  the  sampled  value  is  compared  to  a 
threshold  for  the  purpose  of  making  a detection  decision. 

hi  the  following  sections,  results  of  analyses  relevant  to  detector  per- 
formance are  presented  for  the  three  multiplier  waveforms  of  Figure  3.  Var- 
ious methods  of  comparing  the  sampled  integrator  output  to  a threshold  are 
possible.  The  appendix  presents  some  of  these  for  comparative  purposes. 
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PROCESSOR  ANALYSIS 


Analysis  of  the  qua  si -coherent  video  signal  processor  is  presented 
in  the  following  paragraphs  for  each  of  the  three  multiplier  waveforms.  Equa- 
tions ( 3)  and  (4)  represent  the  signal  input  to  the  I and  Q channel  video 
multipliers.  The  sinusoidal  waveform  of  Figure  3 is  first  considered. 


. \ 1 . 


Figure  3.  Video  multiplier  waveforms:  (a)  sinusoidal  wave,  (b)  square  wave, 
( c)  rectangular  pulse  train. 


A.  Sinusoidal  Video  Multiplication 


The  I and  Q sinusoidal  multiplier  waveforms  are  expressed 
as 

gx(t)  - cos  (0^+  0^  (6) 

and 

gy(t)  - sin  (0)^+  0^  . (7) 

Above,  and  0^  represent  the  best  local  estimate  of  the  doppler  frequency 

and  the  phase,  « and  0.  Subsequent  to  multiplication,  the  signals  in  the  two 
d 

channels  become 

&(t)  =“  c(t)c(t+  t)  {cos  [(a^  - u>d)t  + ^ - #1 

+ cos  [( Wj  + wd)t+  #1}  + c(t+  T)n  (t)  cos  ( to^t  + 0^  (8) 

and 

$(t)  = c(t)c(t+  t)  {cos  [(o^  - «d)t+  ei  - 0] 

- cos  [(c^  + «d)t+  0j  + 0]}  + c(t+  r)ny(t)  sin  ( co^t  + 0^)  . (9) 

The  purpose  of  both  an  I and  a Q processor  channel  can  now  be  seen  by  con- 
sideration of  Equations  (8)  and  (9) . The  first  term  in  each  equation  represents 
the  desired  signal,  whereas  the  second  terms  are  a form  of  self-interference. 
The  latter  terms  actually  represent  the  upper  sideband  resulting  from  the 
mixing  (i.  e. , multiplying)  operation  in  each  processor  channel,  and  they  can  be 
removed  by  summing  the  I and  Q channel  signals.  Further,  it  is  recognized 
that  summing  the  two  signals  provides  coherent  addition  of  the  desired  signal, 
while  the  noise  in  the  two  channels  adds  incoherently,  since  the  noise  com- 
ponents are  Independent.  Defining  the  error  in  the  estimation  of  doppler 
frequency  and  the  phase  as 


and 


(ID 


0^  - 0 - A0  , 

the  signals  of  Equations  (8)  and  (9)  are  summed  to  give 
&(t)  + $(t)  =Ac(t)c(t+  t)  cos  (Awt  + A0) 

+ c(t+ r)  (n^t)  cos  (o^t*  0^  + ny(t)  sin  (0^+ 0^}  . (12) 


Digressing  briefly,  the  Integrator  of  Figure  2 Is  considered  for  the  pur- 
pose of  presenting  its  frequency  response  characteristics.  It  is  noted  that  this 
device  integrates  only  over  a finite  time  Interval;  or  it  performs  the  operation 
of  integrate  and  dump.  With  a gain  factor  1/T,  and  an  Integration  interval  of 
T,  the  integrator  has  the  frequency  response 


I H(  f)  | = 


sin  (7rfT) 
7rf  T 


(13) 


This  function  is  sketched  in  Figure  4.  It  can  be  shown  that  the  two-sided  noise 
equivalent  bandwidth  is  1/  T.  Further,  with  white  input  noise  of  two-sided 
power  spectral  density  tj/ 2,  the  output  power  becomes 


N 


0 


_ JL 

2T 


(14) 


Recalling  the  foregoing  discussion  concerning  the  presence  of  self- 
interference in  the  I and  Q processor  channels,  the  following  observation  can 
be  made.  With  the  aid  of  Figure  4,  it  is  seen  that  this  interference  would 
appear  at  the  output  of  the  integrator,  except  in  a case  where  the  frequency 

co  + oj  has  a period  which  is  some  integer  multiple  of  T,  the  integration 
d 1 

interval. 


Returning  now  to  the  input  of  the  integrator  shown  in  Figure  2,  and  to 
Equation  (12),  attention  is  restricted  to  consideration  of  the  in-range  channels. 
The  relation  of  Equation  ( 5)  is  assumed,  and  this  leads  to  the  following 
expression  for  the  integrator  output  at  time  T. 


9 


m 


Z(T)  * T~Z  (8in  (Ao>T+  A0)  - sin  (A0)} 
AuT 


- l 

+ ^ / c(t+  r)  {nx(t)  cos  ( Wjt  + 0X)  + n (t)  sin  (o^t*-  0^} 


» U(A«T,A0) +n  (T) + n (T)  . (15) 

x y 

The  desired  signal  component,  U(Ao>T,  A0),  Is  first  considered.  In  the 
absence  of  frequency  and  phase  estimation  errors, 

lim  {u(Ao)T,A0)}  = A (16) 

AwT— 0 

A0—O 

To  gain  some  Insight  into  the  effects  of  estimation  errors  upon  signal  degrada- 
tion, it  is  assumed  that  both  errors  are  independent  random  variables  uniformly 
distributed  with  probability  densities 


-p  s AuT  - 0 


and 


f( A0)  = - — , -as  A0<a 

2a 


Degradation  in  signal  level  is  now  defined  by 


(17) 


(18) 


U(Aa>T,A0) 

A 

Averaging  on  both  Ao>T  and  A 0 gives  the  average  degradation, 

pa 

D = — ■ f / U(AwT,A0)d(A0)d(AwT) 

-p  -a 

_ sin  (a)  Si(p) 
a p 
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(19) 


(20) 


.trpytj 


with  the  sine  Integral  defined  by 


J 


I 


p » i 

L 
: 1 

Cs 

i 


Prior  to  completing  the  discussion  on  the  effects  of  signal  degradation, 
the  noise  components  in  Equation  ( 15)  must  be  considered.  Multiplication  of 
the  noise  by  the  locally  generated  replica  of  the  PN  code  acts  to  spread  the 
noise  spectrum,  at  best;  and  it  leaves  the  noise  unchanged,  at  worst.  Which  of 
these  two  possibilities  actually  occurs  depends  upon  the  bandwidth  at  the  point 
where  multiplication  takes  place.  To  avoid  this  question,  it  is  assumed  that 
after  being  multiplied  by  the  PN  code,  the  noise  has  two-sided  spectral  density 
of  t)/2,  at  the  input  to  the  video  multipliers  in  each  channel.  This  assumption 
can  be  made  without  loss  of  generality  in  the  analyses.  Next,  multiplication 
by  the  unit  amplitude  sinusoidal  waveforms  in  the  video  reduces  the  noise  power 
density  by  a factor  of  two.  (The  signal  suffers  the  same  loss. ) With  reference 
to  Equation  (14) , the  total  noise  variance  at  the  output  of  the  integrator  is  77/ 2T 
after  summing  the  variances  of  the  I and  Q noise  components.  Further,  the 
ouiput  noise  is  gaussian. 

i 

Data  obtained  from  Equation  (20)  appears  in  Figure  5 in  the  form  of  con- 
tours of  average  degradation.  Although  this  development  does  not  guarantee 
that  detector  performance  degradation  obeys  Figure  5,  it  has  been  shown  via 
computer  simulation  that  these  data  are  valid  at  50%  probability  of  detection, 
and  approximately  so  below  this  value  (Figure  6) . At  higher  detection  prob- 
abilities, corresponding  to  relatively  large  values  of  S/  N ratio,  degradation  in 
detector  performance  is  greater  than  Figure  5 would  imply.  This  is  an  inter- 
esting result,  since  Equation  (20)  was  derived  for  an  infinite  S/N.  However, 

Equation  (20)  gives  the  average  value  of  signal  degradation.  If  values  of  AwT 
and  A 9 are  selected  from  the  set  of  permissible  values  corresponding  to  1-dB 
degradation  (for  example,  Ao>T  = 1 rad  and  A 6 = 0. 6 rad) , it  is  found  using 
Equation  (15)  that  the  signal  component  degrades  7.2  dB  for  this  particular  set 
of  values.  Such  large  values  occur  with  very  low  probability  when  the  average 
degradation  is  1 dB.  But,  at  the  higher  S/  N ratios  where  detection  probability 
exceeds  90% , only  a very  few  detection  misses  are  necessary  to  degrade  detec- 
tor performance  significantly.  That  is,  a significantly  larger  signal  level  is 
required  in  order  to  eliminate  these  detection  misses  which  occur  with  low 
probability. 
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Finally,  with  regards  to  the  use  of  Figure  5 at  low  values  of  S/N  ratio, 
corresponding  to  detection  probability  below  approximately  10% , it  is  recognized 
that  the  noise  plays  an  important  part  in  the  number  of  returns  detected.  Fur- 
ther, it  is  expected  that  the  data  of  Figure  5 will  reflect  reasonably  accurate 
detector  performance  on  the  average  since  many  detection  misses  are  encoun- 
tered. Figure  6 permits  comparison  of  the  performance  of  an  optimum  coherent 
detector  with  that  of  the  quasi -coherent  detector  operating  with  1-dB  average 
signal  degradation.  Probability  of  detection  for  the  optimum  coherent  detector 
is  derived  in  the  appendix,  and  appears  there  as  Case  1. 

B.  Square  Wave  Video  Multiplication 

Digital  signal  processing  has  its  obvious  advantages.  If  the 
quasl-coherent  signal  processor  of  Figure  2 were  digitally  implemented  to 
process  the  video,  it  would  be  expedient  to  use  square  wave  multiplier  wave- 
forms with  values  ±1.  Such  video  multiplications  are  tantamount  to  simply 
changing  the  sign  of  the  video  signal  when  it  is  negative;  i.  e. , the  video  signal 
is  rectified,  but  the  noise  Is  not  rectified  (a  step  which  would  adversely  affect 
detector  performance) . The  effect  of  square  wave  video  multiplication  is  now 
considered. 

The  square  wave  multiplier  waveforms  are  shown  in  Figure  3 as  applied 
in  the  I and  Q channels  of  the  processor.  These  waveforms  can  be  expressed 
as  Fourier  series,  with  the  appropriate  phasing  for  I and  Q.  The  series  are 


n-1 

4 / 2 

gJ4)  = Z £ cos  (nw  t+nfl  ) (22) 

' n=l, 3,5,...  “ 1 1 

and 

gy(4)  " “ Z “ sin  (naijt+n^)  . (23) 

n=l y 3| 5) • • • 

Equations  (3)  and  (4)  are  multiplied  by  Equations  (22)  and  (23),  respectively, 
to  give  the  multiplier  outputs.  Upon  summing  the  two  channels,  certain  terms 
cancel  while  others  add.  The  composite  signal  input  to  the  integrator  is 
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Pi 
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Restricting  attention  to  in-range  channels  where  Equation  ( 5)  holds,  the  inte- 
grator oulput  at  time  T is 


‘<T>-T  £ n[(n  - 1)I"+  A^Tfsln(1(,'-1)“lt  A“1T 

re=l,5,9,...  1 

+ (n  - 1)01  + A0)  - sin  [(n  - l)01  + A0]} 


- T ~~Ti TT^ ; (sin  ([(n  + l)u,  - AulT 

+ (n  + 1)0  - A0)  - sin  [(n  + 1)0  - A0J}  + n (T)  + n (T) 

1 x x y 


U(AuT,A0) + n (T) +n  (T) 
x y 


(25) 


The  higher  harmonics  of  u,  which  appear  (i.e. , 4,  8,  12,  . . .th  harmonics) 

1 4 

are  a form  of  self-interference.  The  power  in  each  reduces  as  1/n  as  a 


function  of  the  harmonic  number.  With  reference  to  the  sketch  shown  in 
Figure  4,  it  is  noted  that  location  of  the  harmonics  at  zeros  of  the  integrator 
frequency  response  appears  advantageous.  That  is,  let 
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where  m is  one  of  the  doppler  channels  appearing  In  Figure  1.  ( For  purposes 
of  the  foregoing  discussion,  m » 1).  If  Equation  (26)  Is  satisfied  for  the  4th 
harmonic,  all  higher  harmonics  will  also  be  located  at  zeros  of  the  Integrator 
frequency  response.  Under  such  conditions,  the  self-interference  is  minimized 

2 

and  that  which  remains  Is  a function  of  AojT. 


Under  the  conditions  of  Equation  ( 26) , the  signal  at  the  output  of  the 
integrator  In  the  absence  of  doppler  frequency  and  phase  estimation  error  Is 
given  by 

lim  {U(A«T,  Afl)}  =»  — . (27) 

AwT— 0 * 

A0—O 

The  subject  of  average  signal  degradation  due  to  frequency  and  phase  estimation 
errors  is  of  interest  here  as  in  the  case  of  sinusoidal  video  multiplication. 

Here,  using  the  same  symbol  as  before  to  denote  dogradatton, 

»-sW11  • (-) 

Again,  assuming  uniformly  distributed  probability  densities  for  these  errors. 
Equations  (17)  and  (18)  give  the  assumed  densities.  Following  the  approach 
indicated  by  Equation  (20) , the  average  signal  degradation  is  obtained  after 
considerable  mathematical  manipulation.  The  resulting  expression  is 

5 = 2/3  E J{co«I(n-l)6  ][a[(n-l)«  T+/8J 

n=l  9 5|9)  • • * 

- Si[(n  - 1)  WjT  - p]  ] - sin  [(n  - 1)0^  [cin  [(n  - lJ^T  + 0] 

-cln  [(n-l)WlT-/9]]} 


The  so-called  self-interference  apparently  occurs  as  either  constructive  or 
destructive  interference.  If  these  two  possibilities  are  equally  likely,  then 
the  restricting  condition  of  Equation  (26)  might  not  be  necessary.  This  point 
has  not  been  investigated. 
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The  sine  Integral  appearing  in  Equation  (29)  is  defined  by  Equation  (21),  while 

. v 3 

cln  ( • ) is  defined  by 

cin  ({)  = / dt  . (30) 

0 1 

Equation  (29)  has  been  studied  under  the  restriction  given  by  Equation 
(26).  Under  such  conditions,  the  average  signal  degradation  is  independent  of 
the  doppler  frequency  while  dependence  upon  the  initial  phase  is  very  small. 

It  is,  therefore,  the  case  that  Figure  5 provides  a very  good  approximation  of 
average  signal  degradation  when  the  conditions  of  Equation  ( 26)  are  satisfied. 

Returning  to  further  consider  Equation  (25),  the  noise  components  there 
are  easily  analyzed.  Since  the  multiplier  waveform  has  equally  probable  values 

4 

of  ±1,  the  presence  of  the  multiplier  does  not  affect  the  noise.  Denoting  the 
two-sided  noise  spectral  density  in  each  video  channel  by  ij/2,  the  variance  of 
the  noise  at  the  output  of  the  integrator  is  rj/T. 

At  this  point,  some  insight  into  the  performance  of  the  detector  using 
square  wave  video  multiplication  can  be  provided,  hi  the  absence  of  degrada- 
tion, the  signal  given  by  Equation  (27)  is  seen  to  be  4/  ir  times  that  given  by 
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This  is  true  only  if  the  frequency  of  the  square  wave  is  less  than  twice  the 
preintegrator  video  bandwidth.  Otherwise,  the  noise  spectrum  is  spread  by 
square  wave  multiplication. 
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Equation  (16),  the  latter  actually  representing  optimum  coherent  detection. 
However,  the  noise  variance  in  the  square  wave  case  is  twice  that  for  coherent 
detection.  These  numerical  values  show  that  the  S/  N ratio  for  square  wave 
video  multiplication  is  0. 9 dB  below  the  same  ratio  obtained  for  the  coherent 
detector,  hi  the  absence  of  estimation  errors,  computer  simulation  of  the 
processor  using  square  wave  video  multiplication  confirms  that  0. 9-dB  degrada- 
tion in  detection  probability  exists.  Simulation  data  were  also  obtained  with 
the  processor  operating  in  the  presence  of  frequency  and  phase  estimation 
errors.  Representative  results  are  presented  in  Figure  7 corresponding  to 
1-dB  average  signal  degradation  (A®  = 0.6  rad  and  Ao>T  = 1 rad) . As  pre- 
dicted, at  50%  probability  of  detection,  performance  is  1. 9-dB  Inferior  to 
optimum  coherent  detection.  Pronounced  deterioration  is  noted  at  the  higher 
values  of  detection  probability. 

Next,  rectangular  pulse  video  multiplication  is  considered  and  compared 
with  the  two  cases  previously  presented.  Due  to  the  similarity  in  the  waveforms, 
the  material  previously  presented  can  be  used  in  presenting  the  next  case.  This 
is  fortunate  for  all  concerned. 


C.  Rectangular  Pulse  Video  Multiplication 

The  rectangular  pulse  waveforms  appearing  in  Figure  3 are 
now  considered  as  the  video  multiplier  waveforms.  The  significant  charac- 
teristics of  this  waveform  Include  the  fact  that  it  (i.  e. , the  notch)  removes 
noise  from  the  input  of  the  integrator  during  portions  of  the  cycle  of  the  video 
signal  where  the  signal  level  is  low  (i.  e. , in  the  vicinity  of  zero  crossings  in 
Figure  8) . In  addition,  the  notch  occurs  such  that  it  aids  in  reducing  the 
effects  of  doppler  frequency  and  phase  estimation  errors.  The  question  to  be 
addressed  first  concerns  the  width  of  the  notch. 

With  reference  to  Figure  8,  where  the  notch  width  is  defined  as  20  rad, 
a first  attempt  to  determine  the  notch  width  is  to  maximize  the  integrator  output 
S/  N ratio  in  the  absence  of  estimation  errors.  The  quantity  to  be  maximized 
is 
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Figure  8.  Video  waveforms. 
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where  a is  the  noise  variance  out  of  the  integrator  when  <f>  equals  zero,  and  the 
factor  in  parenthesis  in  the  denominator. 


R = 


(32) 


is  a noise  power  reduction  factor.  Taking  the  derivative  of  Equation  ( 31)  with 
respect  to  0 and  setting  the  result  equal  to  zero  yields  the  relation 


(33) 
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It  is  easily  verified  that  the  solution  of  this  transcendental  equation  is 


0.405  rad  . (34) 

Following  the  approach  established  for  the  square  wave  case,  the  Fourier 
series  representation  for  the  rectangular  pulse  train  (having  values  of  0 and 
±1)  is  found  to  differ  from  Equations  (22)  and  (23)  by  the  multiplicative  factor 

C = cos  (n0)  , n = 1,3,5,  ...  (35) 

n 

where  n is  the  harmonic  number.  This  factor  carries  through  the  entire  analy- 
sis of  the  previous  case.  That  is,  it  multiplies  the  signal  component  terms  in 
Equations  (24,  25,  27,  and  29).  Therefore,  the  restrictions  of  Equation  (26) 
are  equally  valid,  and  average  signal  degradation  is  again  given  by  Figure  5 
to  a good  approximation. 

For  the  present  case,  the  signal  at  the  output  of  the  integrator  in  the 
absence  of  estimation  errors  follows  from  Equation  (27)  to  yield  (n  = l) 


lim  {U(AwT,A0)}  - — cos(*)  . (36) 

AwT— -0  * 

A0— 0 

Comparing  the  present  case  with  optimum  coherent  detection,  this  scheme  is 
seen  to  be  inferior  by  0. 35  dB  when  using  Equation  (34)  and  the  noise  reduction 
factor  of  Equation  (32)  in  conjunction  with  Equation  (36) . Again,  computer 
simulation  confirms  0. 35-dB  degradation  in  detection  probability  when  the 
estimation  errors  do  not  exist.  Introducing  frequency  and  phase  estimation 
errors,  the  results  of  computer  simulation  show  that  Figure  5 is  useful  in 
predicting  performance  at  50%  probability  of  detection  and  below.  Illustrative 
data  are  presented  by  Figure  9.  Comparison  of  Figures  6 and  9 reveals 
approximately  0. 35-dB  difference  in  performance  also  exists  between  sinusoidal 
and  rectangular  pulse  video  multiplication  when  estimation  errors  are  present. 


IV.  DESIGN  CONSIDERATIONS 

Performance  characteristics  of  the  quasi-coherent  detector  were 
presented  and  assessed  in  Section  m.  Here  design  related  considerations 
arising  from  the  analytical  treatment  are  presented.  In  view  of  the  performance 
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characteristics  of  the  processor  which  employs  rectangular  pulse  video  multi- 
plication, it  is  the  recommended  configuration  of  those  studied.  With  this 
waveform,  digital  techniques  can  be  employed  for  all  video  processing  opera- 
tions, each  of  which  is  simple  and  straightforward. 

At  the  output  of  the  I and  Q channel  filters  in  Figure  2,  A/D  conversion 
is  assumed.  The  required  operation  involving  the  selected  waveform  is  to 
complement  the  digital  bit  stream  (i.  e. , multiply  by  minus  one)  during  time 
intervals  when  the  pulses  are  negative.  Next  the  I and  Q channels  are  summed. 
The  integrator  shown  in  Figure  2 now  becomes  an  accumulator.  Periodically, 
the  accumulator  is  sampled  and  dumped  (i.  e. , set  to  zero) , and  the  sample  is 
compared  with  a threshold  level. 

The  separation  between  doppler  frequency  channels  and  phase  channels 
can  be  determined  through  the  use  of  Figure  5.  In  addition,  a restriction  on 
the  'location'’  of  the  center  of  the  frequency  channels  (i.  e. , the  values,  o>m) 

is  given  by  Equation  (26) . There  is  no  such  restriction  on  phase  if  Equation  (26) 
is  satisfied.  Consideration  of  Figure  5 shows  that  detector  performance  is  less 
sensitivie  to  doppler  frequency  estimation  error  than  to  phase  estimation  error. 
However,  frequency  estimation  error  is  directly  proportional  to  the  integration 
time,  T,  and  it  accumulates  rapidly  for  higher  doppler  frequencies  (e.g. , 

25  kHz) . Several  tradeoffs  are  available.  Frequency  and  phase  channel  sepa- 
ration can  be  varied  along  the  contours  of  constant  degradation  in  average 
signal  shown  by  Figure  5.  There  also  exists  a trade  between  predetection 
(coherent)  integration  and  post-detection  (noncoherent)  integration,  providing 
some  flexibility  in  selecting  the  integration  time,  T,  and  channel  separation. 

This  point  deserves  further  study  in  an  attempt  to  optimize  performance  and 
minimize  the  required  number  of  channels. 

Finally,  it  is  recognized  that  phase  channels  which  collectively  span  2 ir 
rad  provide  a positive  and  a negative  value  of  the  sampled  integrator  output. 

If  8 small  loss  (i.  e. , several  tenths  of  a decibel)  is  acceptable,  threshold  com- 
parison can  be  made  with  the  magnitude  of  the  output  sample.  Such  an  approach 
then  eliminates  half  of  the  phase  channels.  The  remaining  phase  channels 
cover  the  range  O-x  rad. 

V.  SUMMARY  AND  RECOMMENDATIONS 

This  report  presents  the  evolution  of  a video  qua  si -coherent  detec- 
tor having  potential  application  to  quiet  radar.  Results  of  particular  interest 
are  analyses  and  data  pertinent  to  the  comparative  performance  of  the  processing 


schemes  which  employ  sinusoidal  and  rectangular  pulse  train  video  multiplica- 
tion. It  is  shown  that  the  latter  processor  structure  performs  approximately 
0. 35  dB  below  the  former  in  the  presence  and  absence  of  doppler  frequency 
and  phase  estimation  errors.  The  performance  of  both  processors  is  compared 
to  that  of  an  optimum  coherent  detector,  which  provides  a performance  bench- 
mark. Design  criteria  are  provided  for  separating  the  doppler  frequency  and 
the  phase  channels.  The  rectangular  pulse  processor  configuration  is  recom- 
mended because  it  is  conducive  to  the  use  of  digital  techniques  which  offer 
fast,  compact,  and  relatively  inexpensive  computational  capabilities. 

Further  study  of  the  recommended  processor  configuration  is  necessary 
in  several  areas.  These  areas  include  investigation  of  the  operation  under 
out-of-range  conditions,  and  determination  of  the  tradeoffs  available  between 
predetection  (coherent)  integration,  and  post-detection  (noncoherent)  inte- 
gration. Simulation  programs  developed  during  this  reported  work  can  be  used 
to  facilitate  such  future  effort. 


Appendix.  DETECTOR  PERFORMANCE  FOR  VARIOUS  SIGNAL  PROCESSOR 
CONFIGURATIONS  AND  THRESHOLD  COMPARISON  SCHEMES 


Various  video  signal  processing  configurations  and  methods  of  making 
detection  decisions  by  threshold  comparison  are  considered  in  this  appendix 
with  the  detector  performance  expressed  as  the  probability  of  detection  in  each 
case.  Curves  and  data  are  presented  which  depict  detection  probability  as  a 
function  of  S/N  ratio  for  several  values  of  detector  false  alarm  probability. 
These  results  correspond  to  a single  look,  as  post  detection  integration  is  not 
considered  here. 

For  purposes  of  general  usage  in  this  appendix,  the  following  expressions 
of  signal  plus  noise  are  of  interest: 

x(T)  - A cos  ( 0 ) + n (T)  - u(T)  + n (T)  (A-l) 

y(T)  = A sin  (0)  + ny(T)  - v(T)  + ny(T)  (A-2) 

where  T is  the  time  at  which  the  integrator  is  sampled  ( Figure  A-l) . In 
Equations  (A-l)  and  (A-2),  A is  a positive  real  constant,  unless  otherwise 
specified,  and  it  is  representative  of  the  presence  of  a signal.  The  phase,  0, 
is  a random  variable  which  is  uniformly  distributed  in  the  Interval  -jr  to  ir  rad. 
The  terms  n^(T)  and  ny(T)  represent  noise,  and  these  are  assumed  to  be 

independent  gaussian  random  variables  with  zero  mean  value  and  variance, 

2 

cr  , taken  from  narrowband  gaussian  noise  processes.  Further,  the  constant 
A,  the  phase,  and  the  noise  components  are  assumed  to  be  mutually  independent. 
The  probability  density  functions  for  the  noise  and  signal  components  are, 
respectively,  given  by 


f(n  ) = — “ — exp  -n2/2<r2  (A-3) 

X <r  nT 2?  X 


2 COS  <wt) 


2 SIN  tut) 


Figure  A-l.  I and  Q channel  video  signal  processor. 


The  expressions  of  Equations  (A-l)  and  (A-2)  represent  sampled  outputs 
of  a signal  processor  such  as  shown  in  Figure  A-l.  These  experssions  are 
modified  as  necessary  to  fit  the  various  cases  considered  In  the  following 
developments. 

1.  Case  Is  A Coherent  Signal  Processor  and  Optimum  Detection 


hi  evaluating  the  performance  of  a detector,  it  is  useful  to  compare 
its  performance  with  that  of  the  optimum  detector  under  idealized  conditions. 

For  this  reason,  this  case  entails  use  of  a coherent  signal  processor  which 
requires  that  the  received  frequency  and  phase  be  exactly  known.  The  sampled 
processor  output  then  corresponds  to  Equation  (A-l)  with  6 set  equal  to  zero. 

In  Figure  A-l,  this  means  that  the  signal  is  entirely  in  the  I processor  channel. 
Detection  takes  the  form  of  comparison  of  the  sampled  outputs,  x(T),  with  some 
threshold  level  which  produces  a desired  probability  of  false  alarm. 

The  probability  density  function  for  x(T)  is 

f(x)  = 1 - exp  [-(x  - A)2/2(t  2]  , -®  x ® (A-5) 

<7  27r’ 

. 
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as  shown  In  Figure  A-2.  To  establish  the  threshold  level,  tfa,  from  a desired 

false  alarm  probability,  P , the  required  relationship  is 

w 


Values  of  threshold  level  can  be  obtained  from  Figure  A-3,  where  the  false 
alarm  probability  expression  of  Equation  ( A-6)  is  plotted  as  a function  of  the 
normalized  threshold,  t / <r  . 
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With  signal  present,  the  probability  of  detection  Is  expressed  by 


Vp(x>th) " f ^x)dx 

Si 


For  A > t^»  it  Is  noted  that 
erfc(-t)  = 2 - erfc(t) 


(A-8) 


t > 0 


(A-9) 


The  detection  probability  given  by  Equation  (A-8)  Is  shown  In  Figure  A-4  as  a 
function  of 


SNR  = 


(A-10) 


2<r 


/ V 2 

for  several  values  of  false  alarm  probability.  In  Equation  (A-10),  or  is  the 
variance  of  the  noise  out  of  each  integrator  in  Figure  A-l.  For  the  general 
case,  where  0 is  not  zero,  the  SNR  previously  defined  represents  the  S/  N ratio 
following  each  integrator  of  Figure  A-l,  averaged  over  all  values  of  6.  It  also 
represents  the  S/  N ratio  of  the  coherent  sum  of  the  two  channels,  averaged 
over  all  0.  In  any  case,  this  ratio  is  dimensionally  proper  to  define  It  as  an 
S/N  ratio. 

2.  Case  2:  A Coherent  Signal  Processor  with  Magnitude  Threshold 
Comparison 

The  signal  processor  in  this  case  is  that  considered  In  the  previous 
case.  However,  here  detection  decisions  are  based  upon  comparison  of  the 
magnitude  of  the  sampled  outputs  of  the  processor  with  a threshold  level.  The 
probability  density  function  of  Interest  is  that  for 

z(T)  = |x(T)  | , (A-ll) 
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In  the  absence  of  signal,  the  threshold  level  is  related  to  the  false  alarm 
probability  by 


PF-P(|xl  X^)  = / f)x|(x)dx  = 2 / f(x)  dx  , A =*  0 
*h  *h 


erfc 


(A-13) 


It  is  noted  that  the  probability  density  function,  f(x) , in  the  absence  of  signal 
is  given  by  Equation  ( A-3) . The  previously  mentioned  false  alarm  probability 
is  graphically  shown  in  Figure  A-3. 

When  a signal  is  present,  the  detection  probability  becomes 


* iii  > v > / f|xl 

n 


(x)  dx 


+ erfc 


(A-14) 


The  detection  probability  of  Equation  (A-14)  is  plotted  in  Figure  A -4  where  the 
performance  of  the  detection  schemes  of  Case  1 and  Case  2 can  be  compared. 
As  noted,  the  two  methods  compare  favorably  for  large  SNR.  This  is  due  to 
the  fact  that  the  second  term  above  vanishes  under  large  signal  conditions 
since 

lim  erfc(t)  = 0 , 

t— eo 

and  when  A is  much  larger  than  both  threshold  levels,  Equation  (A-14)  reduces 
to  Equation  (A-8).  For  small  values  of  SNR,  the  difference  in  detector  per- 
formance is  seen  to  be  significant. 
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3.  Case  3:  A Modified  Coherent  Signal  Processor  with  Positive  and 
Negative  Threshold  Comparison 

hi  this  case,  the  coherent  processor  of  Case  1 is  assumed  to  be 
modified  such  that  it  produces  a 180*  phase  ambiguity  in  its  estimation  of  the 
received  signal  phase.  Thus,  0 in  Equation  (A-l)  is  0 or  * rad  with  equal 
probability.  Alternatively,  9 can  be  set  to  zero  while  allowing  A to  take  on 
positive  and  negative  values  with  equal  probability.  Here,  two  thresholds, 

±t^,  are  employed  which  are  symmetrical  about  zero.  The  probability  density 
5 

function  is 


f(x)  — {exp  [-(x  - A)2/2q-2]  + exp  [-(x  + A)2/2tr  2]}  . (A-15) 

cr2  V 2ir' 

hi  the  absence  of  signal,  the  relationship  between  threshold  level  and 
false  alarm  probability  is  given  by 

PF  = P(x  > t^)  + P(x  < -t^)  = 2P(x  > t^) 


2 

<jsT2p 


r 2 2 

J exp  (-x  / 2<t  ) dx  - erfc 


(A-16) 


The  folse  alarm  probability  is  graphically  shown  in  Figure  A-3,  as  it  is  the 
same  function  obtained  previously  in  Equation  (A-13). 

The  probability  of  detection  in  this  case  is  given  by  the  probabilities 

PD  = P(x  >th)  + P(x  < -t^) 


5 . 

The  density  function  of  Equation  (A-15)  is  obtained  by  convolving  the  prob- 
ability density  function  for  the  gausslan  random  variable  with  that  for  A.  The 
latter  density  function  is 

fA( «)  = | 6(  u - A)  + | 6(  u + A) 
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Due  to  the  symmetry  displayed  by  the  probability  density  function,  the  two 
probabilities  in  the  previously  mentioned  expression  are  equal.  The  detection 
probability  can  be  written  as 

* "*h 

pDm  J f(x)  dx+  f f(x)  dx - 2 / f(x)  dx 

*h  - 


'2  <erfc 


&7)+ericft?7)}  • 


(A-17) 


with  f(x)  expressed  by  Equation  (A-15).  Equation  (A-17)  is  the  same  exp  res  - 
sion  obtained  for  detection  probability  in  Case  2,  and  given  by  Equation  (A-14). 
Since  the  threshold  levels  are  also  related  to  false  alarm  by  Identical  expres- 
sions, the  two  detectors  perform  in  an  equivalent  manner. 


4. 


Case  4:  A Modified  Coherent  Signal  Processor  with  Magnitude 
Threshold  Comparison 


Here,  as  in  the  previous  case,  it  is  assumed  that  the  coherent 
processor  is  modified  such  that  a 180*  phase  ambiguity  in  the  estimation  of  the 
received  signal  phase  results.  The  detection  decision  is  based  upon  comparison 
of  the  magnitude  of  the  sampled  output  with  a threshold  level.  In  the  absence 
of  signal,  the  situation  is  the  same  as  that  in  Case  2.  Thus,  the  probability 
of  false  alarm  and  threshold  relationship  is  given  by  Equation  ( A-13) . 

The  density  function  for  x(T)  is  given  by  Equation  (A-15) , with  signal 

present.  Accounting  for  the  magnitude  of  x(  T) , the  probability  density  function 
becomes 


'i*iw 


<r  n/  2*' 

0 s x ^ ® 


{exp  l-(x- A)2/2<r2]  + exp  [-(x  + A)2/2o  2]} 


(A-18) 
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The  detection  probability  is,  therefore. 


GO 


(A-19) 


This  result  too  Is  the  same  as  that  obtained  In  Case  2 for  the  probability  of 
detection.  Thus,  the  two  detection  schemes  are  theoretically  equivalent,  and 
the  detection  probability  is  shown  in  Figure  A-4  as  a function  of  SNR. 


5.  Case  5:  A Noncoherent  Signal  Processor  with  Magnitude  Threshold 
Comparison 


A noncoherent  signal  processor  consisting  of  only  the  I channel  in 
Figure  A-l  is  considered  where  the  received  signal  has  unknown  phase  in  the 
presence  of  additive  gaussian  noise.  The  sampled  output  of  the  processor  is 
described  by  Equation  (A-l) , where  the  parameter  A is  a positive  real  con- 
stant, 9 represents  a random  variable  uniformly  distributed  in  the  range  -it 
to  7 r rad  and  n ( T)  represents  a zero  mean  gaussian  random  variable  with 
2 ^ 

variance,  a . With  reference  to  Equation  (A-l) , the  probability  densities  for 
nx(T)  and  u(T)  are  given  by  Equations  (A-3)  and  (A-4),  respectively. 

Assuming  independent  phase  and  noise,  the  probability  density  function  for 
x(T)  can  be  expressed  in  terms  of  the  characteristic  functions  for  the  two 
densities. 


Concerning  f(u) , the  characteristic  function  is 

M ({)-  / f(u)  exp  ( u)  du  = ~ / ?ffi”g«du 

-a>  -A  N A - u 


Making  the  variable  change,  u = A cos  {<f> ) , the  integral  can  be  evaluated  to 
produce 


where  J ( * ) Is  the  zero  order  Bessel  function  of  the  first  kind.  Similarly!  the 
0 

characteristic  function  for  the  gaussian  random  variable  with  probability  density 
f(n)  is 

00 

M (|)  = — - — • f exp  (-n2/2 a ) exp  (j|n)  dn 

o nT2F  -® 

= exp  (-|2or2/2)  . (A-21) 

hi  terms  of  the  two  characteristic  functions,  the  density  function  for  x(T)  is 
given  by 

00 

f(x)=~  / J0(A|)expU  o2/2)expHSx)d*  , (A-22) 


which  is  identified  with  the  Fourier  transform  of  the  product  of  the  character- 
istic functions.  This  Integral  cannot  be  evaluated  in  closed  form.  One  method 
of  obtaining  numerical  data  representing  this  probability  density  is  to  utilize 

the  fast  Fourier  transform  algorithm  (FFT).6 

The  detector  addressed  makes  decisions  based  upon  the  magnitude  of  the 
processor  output  samples.  Since  both  characteristic  functions  are  even  func- 
tions, the  density  function  Equation  (A-22)  is  real  and  even.  Thus,  the  prob- 
ability density  for  the  magnitude  of  x(  T)  can  be  expressed  in  terms  of  Equation 
(A-22)  as 


Note  that  Equation  (A-22)  must  be  expressed  in  the  proper  form  when  using 
the  FFT  to  evaluate  the  integral.  With  the  discrete  Fourier  transform  (DFT) 
defined  by 

N-l 

G(k)=—  £ g(m)  exp  (-j2nkm/N)  , 

m=0 

Equation  (A-22)  must  be  rewritten  using  the  variable  change,  0 = / 2ir.  Here, 

T is  the  range  of  | over  which  the  characteristic  functions  exist,  and  N is  the 
number  of  sample  points  entered  in  an  N-point  FFT. 


where  the  probability  density  corresponds  to  the  absence  of  signal.  Under  this 
condition,  the  false  alarm  probability  is  again 


[ 


' 

* 

v 


p m — - f exp  (-x2/2 a 2)  dx  ■ erfc  ( • (A-24) 

F J W 27/ 

n 

This  probability  appears  graphically  in  Figure  A-3  as  a function  of  the  normal- 
ized threshold,  t.l <r . 

n 

Li  the  presence  of  signal, 

Equation  (A-23)  to  yield 

CO 

Pp  " P(  lx  I > tfa)  “ 2 / 

*h 

Equation  (A-25)  has  been  numerically  evaluated  using  the  FFT  computations 
of  Equation  (A-22)  to  produce  data  from  which  the  curves  of  Figure  A-7  were 
prepared  for  several  values  of  false  alarm  probability. 

6.  Case  6:  A Noncoherent  I-Q  Signal  Processor  with  Square  Root 
of  the  Sum  of  the  Squares  Threshold  Comparison 

A noncoherent  signal  processor  consisting  of  in-phase  and  quadrature 
channels  as  illustrated  in  Figure  A— 1 is  considered.  The  received  signal  has 
unknown  phase  in  the  presence  of  gaussian  noise,  and  the  sampled  outputs  from 
the  I and  Q channels  are  expressed  by  Equations  (A-l)  and  (A-2) , respectively. 
Further,  the  detection  decision  is  based  upon  comparison  of  the  square  root  of 
the  sum  of  the  squares  of  the  sampled  outputs  of  the  I and  Q channels  with  a 
threshold  level.  Therefore, 

r(T)  = [x2(T)  + y2(T)]1/2  (A"26) 

and  the  probability  density  function  of  interest  is  denoted  by  f(r) . 

It  is  noted  that  the  random  variables  which  represent  the  signal  com- 
ponents of  Equations  (A-l)  and  (A-2),  u(T)  and  v(T),  are  not  independent 
since  they  are  related  by  0.  However,  it  is  important  to  recognize  that  the 
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the  probability  of  detection  is  obtained  from 
f(x)  dx  , 0 ^ x « ® . (A-25) 


density  function  for  Equation  (A-26)  is  independent  of  the  values  assigned  to  0. 
That  is,  the  noise  component  can  always  be  defined  to  be  in-phase  and  in  quad- 

7 

rature  relative  to  any  phasor  representation  of  the  signal  component.  Thus, 
any  convenient  constant  value  of  6 can  be  assumed  which  will  facilitate  deter- 
mination of  f(  r) . Let  0 equal  zero,  and  the  joint  probability  density  function 
for  the  I and  Q channel  outputs  becomes 

f(x,y)  * f(x)f(y)  = 1 — exp  [-(x  - A)2/ 2a  2]  exp  (-y2/ 2 a 2) 

2iro 


Using  the  transformation  of  random  variables 


r 


, 2 2. 1/2 

(x  +y  ) 


x = r cos  ) 
y * r sin  (<f> ) 


r>  0 

-ir  ^ <p  ^ ir 


and 


dx  dy  = r dr  cty 

leads  to  the  relation 

f(r,£)=rf  ( r cos  <f> , r sin  0 ) 
xy 

The  desired  density  function  is  given  by 


7T 

f(r)  = / f(r,0)  d<> 

-7T 

7T 

= — ~ f exp{-[(r  cos«  -A)  +(rsin«)  j/2a  } &f> 
2 iro  -ir 
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Beckmann, 


P.,  Probability  in  Communication  Engineering,  p.  122. 


. Taes  : 
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~ l-(r2  + A2)/2o- Z]  f exp  [ ( rA  cos  * ) /or Z]  d* 


2mr 


-IT 


~expl-(r2  + A2)/2<r2]  IQ 


(A-27) 


where  IQ(- ) Is  the  zero  order  modified  Bessel  function.  The  Equation  (A-27) 

is  often  referred  to  as  the  Rice-Nakagami  density.  As  a matter  of  formality, 
the  joint  density  function,  f(x,y)  could  be  defined  as  a function  of  9 also. 

This  would  lead  to  a need  to  average  over  all  9 in  order  to  obtain  f(  r) . But 
with  9 assumed  uniformly  distributed  in  the  range  -ir  to  x rad,  the  resultant 
would  be  Equation  (A-27).  The  density  function  given  by  Equation  (A-27)  is 
graphically  illustrated  in  Figure  A-8. 


The  false  alarm  probability  from  which  a threshold  level  is  established 
is  obtained  from  Equation  (A-27)  by  setting  A to  zero.  Under  this  condition, 
Iq(0)  has  unit  value,  and  the  Rayleigh  density  function  results.  Thus, 


P(r>tfa) 


; 

-•h 


exp 


$ 


f(r)  dr 


A = 0 


(A-28) 


This  function  is  shown  in  Figure  A-3,  from  which  the  normalized  threshold 
level  can  be  found  for  a given  false  alarm  probability. 

hi  the  presence  of  signal,  the  probability  of  detection  is  expressed  as 


pD-  p(r>  tfa)  - / f(r)  dr 
*h 


(A-29) 


This  Integral  was  numerically  evaluated  to  obtain  data  from  which  the  curves 
of  Figure  A-7  were  constructed. 


•~r  . ■* 
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Figure  A-8.  Probability  density  function,  Case  6. 

In  the  foregoing  development,  it  was  asserted  that  the  received  signal 
phase  could  be  assigned  any  value  for  purposes  of  finding  the  density  function, 
f(r) . This  claim  was  substantiated  during  this  effort  by  digital  computer 
simulation  of  Equation  (A-26),  withx(T)  and  y(T)  given  by  Equations  (A-l) 
and  (A-2).  Several  simulations  were  conducted.  First,  the  phase,  0,  was 
randomly  selected  in  a uniform  manner  over  the  range  of  values,  -7r  to  7r  rad, 
and  the  simulation  was  carried  out.  hi  addition,  0 was  assigned  values  of  10°, 
30*  and  40°,  and  three  separate  simulations  were  carried  out.  Each  simulation 

-4-5-6 

conducted  produced  data  for  false  alarm  probabilities  of  10  ,10  , and  10 

as  a function  of  SNR.  Comparison  of  these  data  with  the  detection  probability 
curves  of  Figure  A-7  showed  excellent  agreement. 
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7.  Case  7 : A Noncoherent  I-Q  Signal  Processor  with  the  of  the 
Magnitudes  Threshold  Comparison 

The  I-Q  signal  processor  of  the  previous  case  is  again  considered 
with  detection  decisions  based  upon  comparison  of  the  sum  of  the  magnitudes 
of  the  sampled  outputs  from  the  I and  Q channels 


z(T)  = |x(T)  | + |y(T)  | 


(A-30) 


with  a threshold  level.  Again,  x(T)  and  y(T)  are  expressed  by  Equations 
(A-l)  and  (A-2) , with  0 uniformly  distributed  In  the  range  -*•  to  ir  rad.  The 
density  functions  for  the  magnitudes  of  x(T)  and  y(T)  are  both  of  the  form  of 
Equation  ( A-22) , but  the  two  random  variables  are  not  Independent. 

In  the  absence  of  signal,  the  probability  density  function  for  z(T)  Is 
obtained  by  convolving  the  magnitude  gausslan  noise  density  with  Itself,  since 
the  noise  in  the  two  channels  is  Independent.  The  density  ftmctlon  for  each  noise 
component  has  the  form  of  Equation  (A-3).  Accounting  for  the  magnitude  of 
each  gausslan  random  variable, 


fU)*  / 


2 / exp  (-n2/2(r2)  exp  [-(x -n)2/2e  2J 

2mr  0 


— ^ exp  (-z2/4<r2)  erf(z/2a)  , 0«z 

tr  nTtT 


Above,  the  error  ftmctlon  Is  defined  by 


(A-31) 


erf(t)=  — J exp  (-r2)  dr 
nTF1  o 


(A-32) 


The  probability  of  false  alarm  Is 

• *h 

PFaP(z>th)a  ft*)**- 1-  / f(z)dz  . 
*h  0 


% 


Substituting  Equation  (A-31)  Into  the  rightmost  integral  of  the  previously 
mentioned  equation,  and  making  the  variable  change,  r ■ z/2<r  , yields 

4 *»/2" 

P = 1 f exp  (-r  ) erf(r)  dr 

F nTT  o 

This  integral  can  be  evaluated  using  the  method  of  integration  by  parts  in  con- 
junction with  the  relation  of  Equation  (A-32)  between  the  two  factors  in  the 
Integrand.  Carrying  out  these  steps  gives  the  false  alarm  probability 

*F-l-^(S  ' <A'33) 

The  function  given  by  Equation  (A-33)  is  shown  in  Figure  A-3  for  purposes  of 
determining  a desired  threshold  normalized  to  a . 

Derivation  of  an  expression  for  the  detection  probability  requires  that  the 
density  function  for  z(T)  be  known.  lacking  such  an  expression.  Equation 
(A-30)  was  simulated  on  the  digital  computer.  The  phase  appearing  in  Equa- 
tions (A-l)  and  (A-2)  was  randomly  selected  from  a uniform  distribution  in 
the  range  of  values,  -it  to  -it  rad.  The  I and  Q noise  samples  were  independently 
selected  from  two  independent  gaussian  distributions.  Results  of  the  simulation 
appear  in  Figure  A-9,  which  depicts  detection  probability.  Each  data  point  in 
the  figure  is  the  result  of  300,000  simulated  samples  of  z(T) . 

8.  A Noncoherent  I-Q  Signal  Processor  with  the  Larger  of  the 
Magnitudes  of  I and  Q Threshold  Comparison 

Once  again,  the  I-Q  channel  processor  of  Case  6 is  considered. 
Here,  detection  decisions  are  made  by  comparing  the  larger  of  the  magnitude 
of  the  I and  Q channel  sampled  outputs 

z(T)  = Max  { |x(T)  I , ly( T)  1}  (A-34) 

with  a threshold  level.  The  probability  density  functions  for  the  magnitudes  of 
x(T)  and  y(T)  are  both  of  the  form  of  Equation  (A-22);  however,  the  random 
variables  for  the  I and  Q samples  are  not  independent. 


In  the  absence  of  signal,  the  I and  Q noise  samples  are  independent. 

8 

Under  this  condition,  the  density  function  is 
fU)-2f|x|<Z)F|jt|U)  . A-0 

with  the  cumulative  probability  distribution 

? 2 ? 2 2 
F,  .(x)  = J f.  ,(x)  dx  = — — J exp  ( -r  /2c  ) dr 

1X1  0 1X1  v'fZiT  0 

= erf  / — - — \ , 0^x5® 

yJTa) 

l 

Thus,  the  density  function  is 


This  function  is  shown  in  Figure  A-3,  from  which  values  of  normalized 
threshold  levels  can  be  obtained. 


As  in  the  previous  case,  the  detection  probability  as  a function  of  SNR 
for  several  values  of  false  alarm  probability  was  found  by  digital  computer 
simulation,  hi  this  case,  the  resulting  data  values  were  found  to  be  in  com- 
plete agreement  with  those  data  obtained  in  Case  7.  Thus,  the  detection 

Q 

schemes  of  Cases  7 and  8 are  equivalent,  a result  stated  by  Nathanson. 
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